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We discuss how the Lagrange multiplier method of non-equilibrium steady state statistical me- 
chanics can be applied to describe the electronic transport in a quantum wire. We describe the 
theoretical scheme using tight-binding model. The Hamiltonian of the wire is extended via a La- 
grange multiplier to "open" the quantum system and to drive current through it. The diagonalization 
of the extended Hamiltonian yields the transport properties of wire. We show that the Lagrange 
multiplier method is equivalent to the Landauer approach within the considered model. 



In recent years, several experimental groups have re- 
ported the measurements of the transport properties of 
carbon nanotubes, self- assembled monolayers of conju- 
gated polymers and even individual molecules 1 . This 
development has attracted much attention from theo- 
reticians and most theoretical works reported so far are 
based on the Landauer type theory 2 . Within the Lan- 
dauer approach it is assumed that incoming electrons 
are scattered along the molecular wire and the con- 
ductance can be directly evaluated from the scatter- 
ing T-matrix. The nonequilibrium Green's functions or 
Lipmann-Schwinger equation approaches have been used 
to compute T-matris^ai^. This extensive calculation 
work on molecular junctions have provided basic under- 
standing of the physical processes although the aim of 
the direct reproducing of experimental I-V characteris- 
tics is still elusive^. This current status of electron trans- 
port calculations in molecular wires requires that other, 
"non-Landauer-type" theoretical descriptions should also 
be properly explored. 

A quantum wire with direct, static current is a 
nonequilibrium steady state system. This can be un- 
derstood based upon the following qualitative discussion. 
Suppose that the left end of the wire works as a equi- 
librium source of the electrons and the right end works 
as a drain. The source and drain have different time- 
independent chemical potentials (fi s > /i^). The source 
pumps electrons into the wire trying to establish the equi- 
librium, but at the same time the drain keeps pulling 
electrons out of it also trying to equilibrate the wire with 
itself. The equilibrium is never achieved as long as the 
chemical potentials fi s and fid are kept fixed or if the 
current is constrained to be time-independent constant. 
The system stays in nonequilibrium steady state provided 
that there is the time-independent current through it. 

Having established that a quantum wire with direct 
current is a steady state nonequilibrium quantum system, 
we can apply the powerful machinery of modern theo- 
retical methods available in nonequilibrium steady state 
statistical mechanics to describe it. Recently, there have 
been considerable advances towards the formal under- 
standing of nonequilibrium steady state open quantum 
systems on the basis of the Lagrange multiplier based 
methoc&Si^iiSi. The basic assumption of the Lagrange 



multiplier based method is that the origin of current is 
irrelevant to physics of the steady state or in other words 
that a homogeneous current-carrying state is the same 
whether it is introduced by the reservoirs, i.e. via bound- 
ary conditions, or by a bulk driving field. A nonequi- 
librium steady state can be produced by the following 
three-step "algorithm"ii&: 

1. Define the operator of current J via the continuity 
equation. 

2. Extend the Hamiltonian by adding the term (— AJ) 
where A is a Lagrange multiplier. 

3. Diagonalize the extended Hamiltonian. 

The aim of this paper is to place electronic transport 
calculations into the general framework of the modern de- 
velopment of nonequilibrium steady state statistical me- 
chanics. We begin with the model tight-binding Hamil- 
tonian for an infinitely long quantum wire 

H = -t 22(al +la a n(T + al ia a n+la ) , (1) 
no- 
where t is the hopping matrix element. The operators 
aJ lCT (aw) create (annihilate) single electron on the site 
n, where a = ±i is the spin of the electrons. Being 
fermionic operators the operators a' na (a na ) obey the 
standard anticommutation relations: 

{^ ncr , G m(J / } — fonm^aa' j (2) 

(4a > a L<r> } = > °w } = ■ ( 3 ) 

The left part of the wire is considered to be source for 
current and the right end serves as drain. We do not 
use any assumptions regarding the physical nature of 
the source/drain we merely assume they exist and are 
able to maintain a steady, i.e. time-independent, current 
through the wire. 

As the first step, we define the operator of current via 
the continuity equation. The number of electrons on the 
site n is given by the expectation value of the operator 

N n = aL-«W ■ (4) 



By making use of Heisenberg representation the continu- 
ity equation can be written as the Heisenberg equation- 
of- motion for the operator N n : 

N n = i [H , N n ] . (5) 

Given the standard anticommutation relations between 
electron creation and annihilation operators, the r.h.s. 
commutator (JjJJ can be readily computed: 

(7 

+^X]( a ^ a "- lCT ~ a n-l<r°W) ■ ( 6 ) 
a 

Comparison of eq. JBJ with the finite difference expression 
for continuity equation N n — ~(j n — jn-i) yields the 
definition of the operator of current through the site n: 

jn — H ^ ] { a n+lo- a n<? ~ 0'n rT a n+la) ■ 
cr 

By making the sum of on-site currents along the wire we 
obtain the net current through the wire 

J = it^iai+l^na ~ a na a n+ l a ) . (7) 
na 

The next step in our scheme is to fix the net current 
via a Lagrange multiplier A. To this end the Hamiltonian 
Hq is modified by adding the term which constraints the 
macroscopic current J: 

H = H — A J = —t'^2( a li+la a na- + a na a n+la) 
n a 

-\it^2(al +la a na - ai a a n+la ) . (8) 

ntj 

The Hamiltonian JSJ is hermitian although it is no longer 
a real operator. The term (—A J) breaks the symme- 
try between electrons moving along the wire in opposite 
directions. Now all physical properties of the system in- 
cluding transport characteristics should be obtained from 
the diagonalization of the Hamiltonian H not Hq . 

The final step is the diagonalization of the extended 
Hamiltonian. Using the periodic boundary conditions 
Q>n<r = o-n+Ncr (N -the "length" of the box) the Hamilto- 
nian (JHJ can be diagonalized via the Fourier transforma- 
tion 

a na = -7= ^ exp(-^fcn)a fcCT , (9) 

^ k 

where the sum over the k runs over the first Brillouin zone 
(k = =£i with i = 0, JV — 1). Furthermore, the oper- 
ators OuAcLka) still obey the fermionic anticommutation 
relations. The transformation JHJ) brings the Hamiltonian 
© to the diagonal form 

H = Y J E{k)al a a ka , (10) 
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k 

FIG. 1: The current-dependent band energies of the one 
dimensional conductor for various current densities: solid line 
- zero current, dashed - 0.025 a.u. and dotted -0.05 a.u. 

with the current-dependent dispersion relation for the 
quasiparticle energy (band energy): 

E(k) = -2i(cos(fc) + Asin(fc)) . (11) 

As we let the current tend to zero, i.e. A — > 0, we recover 
the usual result for the band energy E(k) = — 2tcos(fc). 
The dispersion relation l|llf> is not yet in a form allowing 
for the energy to be computed as the Lagrange multi- 
plier A is not known yet. The additional equation for the 
Lagrangian multiplier A is obtained if the density of the 
expectation value of the net current operator (JJJ over the 
ground state of the Hamiltonian (JSJ 

^ J ) = |E sin ( fc K (12) 

is required to yield the a priori known current density 
/. The occupation numbers n k equals to 1 if the band 
energy E(k) < and zero otherwise for half-occupied 
conductance band, i.e. n% = 8(—E(k)) where 9 is the 
Heaviside step function. 

We finally demonstrate application of the method with 
numerical and analytical examples. The numerical cal- 
culations were carried out for half-occupied conductance 
band with the box length N = 200. The value of the 
hopping integral t is chosen to be 2.5 eV which is the 
standard value for monoatomic metallic wiresii. The La- 
grange multiplier A is obtained via the numerical solution 
of the following nonlinear equation: 

sin(fc)n fe = /, (13) 

k 

In Fig.l, the current-dependent band energies are 
shown for three different values of the current density. 
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Since the commutator [H, J] vanishes exactly it is per- 
fectly possible to characterize each electronic fc-state by 
the value of the current in this state: 

J k = {a ktT Ja\ a ) = 2tsm{k). (14) 

It enables us to distinguish the electronic states by their 
negative or positive contributions to the net current. 
From Fig.l we see that the current re-arranges the band 
structure in such a way that the part of the band which 
corresponds to states with positive current (0 < k < ir) 
has lower energy than the zero-current band while the 
energy of states with negative current are increased with 
respect to the zero-current solution. This selective to cur- 
rent alternation of the band energy splits the occupation 
numbers of states with positive and negative Jk- 

Now we turn our attention towards how the voltage 
drop can be defined within the Lagrange multiplier based 
transport theory. For the most regimes observed in 
molecular wires the electrons are transmitted through the 
levels in the vicinity of the Fermi-energy Ef. Therefore 
to associate voltage with current we need to get clear 
physical picture of Fermi-energy response upon current. 
In the first Brillouin zone there are always two k which 
correspond to the Fermi-energy 

E f = -2f (cos(fe) + A sin(fc)). (15) 

First solution of ea. fl"5|l . k + (0 < k + < n), gives a posi- 
tive contribution to total net current while the second 
solution, fc_ (tt < fc- < 27r), contributes negatively. 
By turning off the current (i.e. in the limit of A = 0) 
the energy of of fc+-state is increased by the value of 
(— 2iAsin(fc + )) while the energy of k- -state goes down 
by (— 2tX sin(fc_)). In the terminology of the Landauer 
theory for reflectionless contact, the energies of states 
originating from the source, i.e. k + states, are effectively 
decreased by voltage while the energies of states occupied 
by electrons from the drain, i.e. fc_ states, are increased 
by voltage. Given that the voltage does the same job as 
the term (— 2iAsin(fc_)) we arrive to the following equa- 
tions: 

E f = -2tcos(fc±) - 2tAsin(fc±)) 

= -2icos(fc±)T^f/, (16) 

where U denotes the voltage drop. Eq.(15) yields the 
following expression for the voltage U : 

U = 2t\(ajx(k+) - sva(k-)) (17) 

The voltage produces splitting ( as in standard Landauer 
theory) the Fermi-Dirac occupation numbers of the the 
states with positive and negative Jk- The Fermi-Dirac 
occupation numbers are given in Fig. 2 for two different 
values of the currents: / = 0.01 a.u. (the lower panel) 
and I = 0.025 a.u. (the upper panel). The Heaviside step 
function is smeared out as the corresponding Fermi-Dirac 
function taken at temperature T = 0.3 eV. The values of 
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FIG. 2: The Fermi-Dirac occupation numbers of current car- 
rying states. The solid line is for the states with positive 
contribution to the net current (Jk > 0). The dashed line 
is for the occupation probabilities of the states with negative 
Jk- Eo(k) — —2tcos(k) is the band energy of the zero-current 
wire. 

the voltage computed by the formula Ijl7(l are also plotted 
on Fig. 2. We would like to emphasize that the definition 
of voltage described above is applicable directly only if 
current J is an integral of motion for the Hamiltonian 
H , i.e. [Hq, J] = 0. The commutator [H ,J] does not 
necessarily vanish exactly in general case and the occu- 
pation numbers of "upstream" and "downstream" elec- 
trons are tangled because the eigenstates of the extended 
Hamiltonian H are no longer uniquely characterized by 
value of the current in these states. Therefore it might be 
practically cumbersome to extract the voltage for a real 
molecular wire by looking at the populations of states 
carrying currents in the different directions. The rigor- 
ous definition of the applied voltage is generally a very 
controversial issue if one does not invoke to noninteract- 
ing electron reservoirs to represent boundary conditions 
and similar problems are encountered in ref^i One can 
always resort to the exact definition of the voltage based 
upon the following physical picture: a voltmeter deter- 
mines the voltage drop along the wire by measuring the 
work required to move a point unit charge from one side 
of the wire to the other. This brute force algorithm of cal- 
culating the work done on a test point charge can be al- 
ways applied although it is not computationally the most 
tractable scheme. The evaluation of the different voltage 
definitions for the correlated electrons will be discussed 
elsewhere. 

To demonstrate a compatibility with the Landauer 
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approach we re-derive the Landauer expression for the 
conductivity 2 " using our formulae for the net current i|12|) 
and our definition of the voltage El Sum (|12|l can be 
replaced by the integral provided that for large N the 
spectrum E(k) can be considered as a quasi-continuum: 



I = 



2t 



2t 



dk sm(k)n,k 



dfcsin(fc)(n fe - n k+n ). 



(18) 



Assuming the small current, i.e. small A, and by mak- 
ing use of the Taylor expansion the following expression 
for the occupation numbers difference can be found 

nk — nk+Tr = 0(cos(fc)) — 0(cos(fc+7r))+2Asin(fc)<S(cos(fe)). 

Only the ^-function term gives a non- vanishing contribu- 
tion to the integral (|18|l and the direct integration results 
into the following simple expression for the current den- 
sity: 



TT 



(19) 



To compute voltage by the formula l|17fl k + and fc_ should 
be determined. For half-occupied band Ef = and the 
ea. (|15fl can be straightforwardly resolved: 



.1 



arctan( — ) 
A 



k + = TT 

Then ea. l|17|) yields the voltage 



fc_ = 2tt- arctan(-) (20) 
A 



u = 



4iA 



Vx 2 



i 



(21) 



Assuming again that the current is small and substituting 
the expressions for current (|19fl and voltage (|21() into the 
definition of conductance G = I/U we obtain: 



G 



1 

TT 



(22) 



which is the Landauer value of G for a single transport 
channel of an ideal one-dimensional lead 12 . Re-deriving 
the standard result we demonstrate that our method is 
equivalent to the Landauer approach within the consid- 
ered model. It should be noted, however, that there is 
no guarantee that this formalism will exactly reproduce 
the Landauer results in more complicated cases. 

Finally, we discuss how the Lagrange multiplier based 
transport theory can be applied to a molecular wire with 
inhomogeneous electron density. To describe this realistic 
scenario we specify the constraint on the current density 
distribution in the following form (constrained continuity 
equation): 



dydzj x {v) = I, 



(23) 



where / is a desired value of the current through the 
wire. Within this description, net current flow is aligned 
along the a;— axis and the net current flow across a cross 
section J dy dz j x (r) is constrained, and this quantity is 
readily available experimentally. It leads to the following 
extended Hamiltonian: 



H = H Q - J dxX(x) 



dydzj x (r), 



(24) 



where A(x) is the pointwise Lagrange multiplier in the 
additional space- and current-dependent term. An in- 
homogeneous interacting electron gas has been consid- 
ered based upon a variational analog of the Lagrange 
multiplier based transport theory and within the den- 
sity functional theory one generally needs to solve the 
set of integral self-consistent equations for the pointwise 
Lagrange multiplier A(x) and the current-carrying Kohn- 
Sham orbitalsi^. 

The final comment regarding the general applicabil- 
ity of the method is in due order. The common char- 
acteristic feature of all nonequilibrium steady state sys- 
tems is the constant flux of a certain physical quantity, 
e.g. heat conduction (energy current), diffusion (parti- 
cle current), electrical conductivity (charge current). We 
have discussed the scheme on an example of tight-binding 
Hamiltonian but it can be extended with a slight mod- 
ification to any open steady state quantum system, e.g. 
the Ising model with the energy flixA Constant current 
Hartree-Fock and Kohn-Sham approximations have also 
been recently formulated within the approachiii^. 

In this paper, we have given the Lagrange multiplier 
based formulation of electronic transport problem. We 
discussed the three-step practical algorithm to produce 
current carrying steady states in quantum wires. First, 
we defined the operator of macroscopic current via con- 
tinuity equation. Next, tight-binding Hamiltonian was 
modified by the Lagrange multiplier term to account for 
the steady current. Then, the current-dependent band- 
energies and the occupation numbers were obtained by 
the diagonalization of the extended Hamiltonian. The 
definition of the applied voltage which is compatible with 
the Landauer description was discussed. A sample cal- 
culation on an albeit simple model system produced ex- 
act agreement with results obtained from the Landauer 
theory. While the one-dimensional tight-binding model 
avoids a number of computational difficulties in using 
this approach, it shows how a Lagrange multiplier can 
be applied to describe electronic transport properties of 
quantum wires. 
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